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基于 改进 高 斯 一 拉 普 拉 斯 算 子 的 噪声 图 像 边缘 检测 方法 
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摘 要 : 针对 现 有 梯度 算 子 在 图 像 边缘 检 测 中 存在 的 对 噪声 比较 敏感 的 问题 ,提出 了 一 种 改进 的 高 斯 - 拉 普 拉 斯 算 子 的 
边缘 检测 方法 。 品 声 图 像 中 的 边缘 检测 是 一 项 关键 任务 ,然而 目前 常用 的 几 种 梯度 算 子 ,包括 已 经 提出 的 高 斯 - 拉 普 拉 
斯 算 子 都 没 能 取得 理想 效果 。 提 出 的 方法 对 传统 的 拉 普 拉 斯 边缘 检测 算 子 作 了 改进 ， 并 与 高 斯 滤波 器 相 结合 。 首 先 ， 
应 用 高 斯 滤波 器 来 平滑 图 像 ， 抑 制 噪声 ; 然后 ， 基 于 拉 普 拉 斯 梯度 边缘 检测 器 进行 边缘 检测 ; 最 后 ， 在 Imagenet 数据 
集中 选取 了 10 幅 图 像 进行 了 实验 , 将 提出 的 高 斯 梯度 边缘 检测 器 与 传统 的 边缘 检测 器 进行 比较 。 评估 结果 显示 , 提出 
的 方法 所 获得 的 峰值 信 骂 比 (PSNR) 高 于 对 比 算 法 ， 而 均 方 误差 (MSE) 更 小 。 实 验 结果 表明 ， 提 出 的 方法 在 实际 应 用 中 
能 够 有 效 提高 噪声 图 像 边缘 检测 的 质量 。 
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中 图 分 类 号 : TP391 doi: 10.3969/j.issn.1001-3695.2018.02.0183 


Noise image edge detection based on improved Gauss-Laplace operator 
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Computer SciencekTechnology, Huazhong University of Science & Technology, Wuhan 430074, China) 


Abstract: Aiming at the problem that the existing gradient operator is sensitive to noise in image edge detection, this paper 
proposed an improved Gaussian-Laplacian edge detection method. Edge detection in noisy images was a key task. However, 
several commonly used gradient operators, including the proposed Gaussian-Laplace operator, had failed to achieve the desired 
results. The proposed method improved the traditional Laplacian edge detection operator and combined it with a Gaussian filter. 
First, it used a Gaussian filter to smooth the image and suppress noise. Then it performed edge detection based on a Laplacian 
gradient edge detector. Finally, it selected 10 images in the imagenet data set. It compared the Gaussian gradient edge detector 
proposed with the traditional edge detector. The evaluation results show that the peak signal-to-noise ratio (PSNR) obtained by 
the proposed method is higher than that of the comparison algorithm, and the mean square error (MSE) is smaller. Experimental 
results show that the proposed method can effectively improve the quality of noise image edge detection in practical applications. 
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动 时 ， 零 交叉 消失 在 其 二 阶 导数 的 尺度 表示 中 。 对 于 2D 信和 号 


J 应 用 ， 也 不 会 随 着 刻度 的 增加 而 产生 零 交叉 。 

图 像 由 具有 不 同 灰 度 级 的 像素 组 成 ， 在 灰 度 级 不 连续 性 检 文献 [6] 提 出 了 一 种 图 像 分 割 模型 ， 其 包括 由 canny 和 归 一 
测 中 ， 尽 管 点 和 线 的 检测 非常 重要 ， 但 边缘 检测 却 是 目前 为 止 ”化 的 切割 特征 向 量 的 边缘 检测 。 在 边缘 检测 之 前 ， 先 采用 基于 
最 为 通用 的 方法 。 图 像 的 边缘 检测 是 利用 图 像 一 阶 导数 的 极 值 ”噪声 类 型 的 中 位 数 、 高 斯 或 Frost 滤波 器 进行 预 处 理 。 文 献 [7] 
(梯度 算 子 ) 或 二 阶 导数 过 零点 〈 拉 普 拉 斯 算 子 ) 信息 来 实现 。 “开发 了 一 种 用 于 医学 图 像 的 边缘 检测 技术 ， 以 基于 强度 梯度 和 


的 07。 边缘 检测 是 确定 和 定位 图 像 中 的 尖锐 不 连续 性 的 过 程 。 ”纹理 梯度 特征 来 追踪 解剖 器 官 的 边界 ， 与 传统 的 有 源 轮 廊 模 型 
在 图 像 分 析 中 ， 边 缘 检 测 起 着 
术 之 一 。 高 斯 滤波 器 在 边缘 检测 器 中 起 着 至 关 
高 斯 滤波 器 被 用 于 一 维 信号 


要 的 作用 ， 它 是 传统 的 分 割 技 。” 相 比 ， 提 出 的 模型 能 够 产生 有 效 的 结果 。 文 献 [8] 在 图 像 上 
要 的 作用 后 。 各 种 基于 梯度 的 边缘 检测 器 。Canny 算 子 能 够 产生 好 的 结 
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平滑 。 当 从 一 个 细 尺 度 到 粗 尺度 移 但 是 在 某 些 情况 下 需要 调整 参数 。 文 献 [9] 证 明了 与 拉 普 拉 


斯 也 
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Sobel 边缘 检测 器 相 比 ,形态 边缘 检测 器 在 存在 噪声 的 情况 下 ， 

对 图 像 产生 有 效 的 结果 。 文 献 [10] 提 出 了 基于 图 像 形态 学 运算 

的 边缘 检测 ， 通 过 取 扩 张 与 侵蚀 图 像 之 间 的 差异 来 确定 边缘 。 

文献 [11] 对 Berkeley 数据 集 上 0 的 各 种 边缘 检测 算 子 进行 了 分 

析 ， 在 图 像 与 边缘 检测 之 前 进行 高 斯 滤波 器 预 处 理 ， 产 生 了 优 
结 


异 的 结果 。 

然而 以 上 提出 的 方法 依然 在 不 同 程度 上 存在 对 噪声 比较 敏 
感 的 问题 。 基 于 此 ， 本 文 提 出 了 一 种 改进 的 高 斯 一 拉 普 拉 斯 算 
子 的 图 像 边缘 方法 。 该 算法 对 传统 的 拉 普 拉 斯 边缘 检测 算 子 作 
了 改进 ,并 与 高 斯 滤波 器 相 结合 。 在 Imagenet 数据 集中 选取 10 
幅 图 像 进行 了 实验 。 与 常用 的 梯度 算 子 相 比 ， 取 得 了 令 人 满意 
的 结果 ， 验 证 了 提出 方法 的 可 靠 性 。 


1 ”常用 的 梯度 算 子 

大 多 数 边缘 检测 算法 都 是 基于 图 像 的 导数 或 梯度 
况 下 ， 图 像 易 被 高 斯 噪声 破坏 或 受到 斑点 噪声 的 影响 。 ,因此 ， 
噪声 图 像 中 的 边缘 检测 就 起 着 至 关 重 要 的 作用 。 
边缘 检测 器 的 目的 是 追踪 图 像 中 期 望 的 感 兴趣 区 域 的 边界 
或 轮廓 。2D 函数 8 ( xz, y ) 的 梯度 如 式 〈1) 所 示 。 
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上 述 矢 量 的 大 小 给 出 ， 如 下 所 示 : 
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式 (4) 决定 。 
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边缘 检测 算法 应 仔细 跟踪 边缘 , 必须 消除 由 于 噪声 像素 引 
起 的 假 边 。 以 下 是 在 图 像 的 边缘 跟踪 中 做 出 的 假设 : 

。 ”边缘 像素 的 梯度 比 噪声 像素 更 大 。 

。 ”边缘 大 小 和 方向 沿 着 边缘 缓慢 变化 。 

以 下 对 这 几 种 常用 的 梯度 算 子 算 子 作出 介绍 。 
1.1 Sobel 算 子 

Sobel 算 子 由 一 对 3x3 的 卷 积 内 核 组 成 。Sobel 算 子 中 的 
内 核对 垂直 和 水 平 运行 的 边缘 产生 最 大 值 。 在 输入 图 像 中 单独 


边缘 强度 


梯度 方向 


应 用 内 核 以 在 每 个 方向 ( G, 和 G, ) 上 产生 梯度 分 量 的 不 同 测 
量 ， 然 后 将 内 核反应 组 合 在 一 起 ， 以 找 出 每 个 点 的 梯度 的 绝对 
幅度 和 该 梯度 的 方向 0。 


式 (5) 给 出 。 


1G |F= arctan 图 (5) 
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梯度 大 小 


Ch | inaz aX (iN\ 会 作 期 刊 
代 文 征 ， 等 


图 1 给 出 了 确定 X( 垂 直 ) 和 》 (水 平 ) 方 向 上 梯度 的 Sobel 
边缘 检测 掩 模 。 图 2 是 Lena 图 像 03] 的 Sobel 算 子 边缘 检测 结 


四 
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1 Sobel 边缘 检测 算 子 的 掩 模 


图 2 原始 Lena 图 像 和 Sobel 算 子 边缘 检测 处 理 后 的 Lena 
像 

1.2 ”Roberts 十 字 梯 度 算 子 

Roberts 十 字 梯 度 算 子 是 主要 用 于 计算 图 像 的 二 维 空间 梯 
， 该 方法 简单 而 快速 。 边 缘 检 测 是 确定 和 定位 图 像 中 的 尖锐 
不 连续 性 的 过 程 。Roberts 十 字 梯 度 算 子 属 于 一 阶 微分 算 子 。 对 
于 噪声 及 光照 不 均匀 等 情况 的 存在 ， 不 能 完整 地 勾画 出 物体 的 
边缘 , 因此 不 能 抑制 噪声 。 Roberts 操作 符 中 的 内 核对 运行 在 45° 
的 边缘 产生 最 大 值 。 图 3 为 Roberts 操作 符 的 掩 模 图 。 图 4 是 
Lena 图 像 的 Roberts 十 字 梯 度 算 子 算 子 边 缘 检测 结果 。 
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3 ”Roberts 边缘 检测 算 子 的 掩 模 
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(b) 
图 4 原始 Lena 图 像 和 Roberts 十 字 梯 度 算 子 边 


缘 检测 处 理 后 


的 Lena 图 像 
1.3 Prewitt 算 子 
Prewitt 算 子 类 似 于 Sobel 算 子 , 用 于 检测 图 像 中 的 导 
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水 平 边缘 。 与 Sobel 操作 符 不 同 ， 
素 。 类 似 于 Roberts 十 字 梯 度 算 子 ，Prewitt 算 子 也 属 卫 


分 算 子 ， 后 


已 不 重视 靠近 撼 模 中 心 的 像 


单 对 噪声 比较 敏感 ， 
Prewitt 梯度 算 子 的 边缘 检测 符 掩 模 如 


噪声 的 影响 。 
图 6 为 Lena 


丽 


了 6 原始 Lena 


传统 的 拉 普 拉 其 


斯 噪声 破坏 或 受到 斑点 


到 像 的 Prewitt 梯度 算 子 边缘 检测 结果 。 


Gy 


于 像 和 Prewitt 算 子 边缘 检测 处 理 后 的 Lena 


基于 改进 高 斯 一 拉 普 拉 斯 算 子 的 边缘 检测 


会 ”2.1 拉 普 拉 斯 算 子 
区 别 于 前 画 
像 的 二 阶 导轨 


了 (x, y) 的 区 


的 几 种 梯度 算 子 ， 拉 普 拉 
来 找到 边缘 并 搜索 过 零点 。 具 有 像素 强度 


银 的 拉 普 拉 斯 变换 工 ( x,y》) 如 式 〈6) 所 示 . 
07 67 
A 
Or Oy 


普 拉 斯 算 子 常用 


边缘 检测 时 ， 拉 普 拉 


产生 双 像素 宽 的 边缘 。 
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图 8 拉 普 拉 斯 算 子 拓 展 模 板 


(b) 


图 9 原始 Lena 图 像 和 拉 普 拉 斯 算 子 边缘 检测 处 理 后 的 Lena 


2.2 高 斯 - 拉 普 拉 斯 算 子 
2.1 节 可 知 ， 拉 普 拉 其 


行 边缘 检测 之 前 ， 先 


| 


于 较 瞳 的 区 域 


更 其 变 得 更 亮 。 


算 子 一 样 ， 拉 普 拉 


了 改进 ， 在 已 知 边缘 


到 像 的 噪声 。 


8 为 改进 的 拉 普 拉 


进 拉 普 拉 斯 算 子 的 边缘 检测 结果 。 


像素 的 暗 区 还 是 亮 区 
斯 算 子 拓展 模板 。 9 为 Lena 攻 


加 


7 拉 普 拉 斯 算 子 常 


模板 


图 像 


i 算 子 并 不 能 抑制 噪声 ， 如 果 有 
种 算 子 能 够 将 高 斯 平滑 滤波 器 和 拉 普 拉 斯 算 子 结合 起 来 ， 在 进 


平滑 掉 噪 声 ， 一 定 能 够 取得 良好 的 效果 。 


甘于 此 ,提出 了 高 斯 一 


斯 算 子 .图 10 所 示 是 常用 的 5x5 


算 子 边缘 检测 输出 。 


的 高 斯 一 拉 普 拉 斯 模板 。 图 11 为 Lena 图 像 的 高 斯 一 拉 普 拉 斯 
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(a) 


10 ”高 斯 - 拉 普 拉 斯 边缘 检测 算 子 模板 


(b) 


11 原始 Lena 图 像 和 高 斯 一 拉 普 拉 斯 算 子 边缘 检测 处 理 后 


2.3 ”高 斯 平滑 滤波 器 


拉 普 拉 斯 边缘 检测 包括 两 


的 Lena 图 像 


级 边缘 增强 和 跟踪 上 9。 应 用 高 斯 


滤波 器 来 平滑 图 像 ,其 较 大 的 


在 确定 边缘 方向 之 后 ， 


核心 尺寸 降低 了 对 噪声 的 敏感 度 。 
最 大 抑制 来 追踪 边缘 的 路 径 ， 并 
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忽略 那些 不 是 边缘 部 分 的 像素 。 最 后 ， 应 用 滞后 姜 值 来 消除 条 
纹 。 定 义 11 > 妇 的 两 个 阔 值 ( 红 和 1 )， 并 应 用 于 图 像 的 梯度 
大 小 。 闵 值 大 于 九 的 像素 被 认为 是 边缘 像素 ， 并 且 连 接 到 大 于 
12 的 边缘 像素 的 像素 被 推定 为 边缘 像素 。 
高 斯 梯度 使 用 高 斯 方法 的 导数 来 确定 标量 2D 图 像 和 3D 
体积 的 梯度 /导数 。 
高 斯 核 的 二 维 表示 如 下 : 


] 过 | 
2 
Nl 二 e~” (7 
高 斯 滤波 器 中 的 O 术语 称 为 平滑 刻度 。 该 比例 对 高 斯 滤波 
器 的 响应 有 很 大 的 影响 。O 的 值 越 大 ， 图 像 将 变 得 模糊 ， 并 且 
对 噪声 的 敏感 性 降低 。 图 12 为 Lena 图 像 的 高 斯 模糊 效果 。 


(b) 
14 ” 灰 度 处 理 的 Lena 图 像 和 高 斯 滤波 处 理 的 Lena 图 像 
b) 高 斯 核 治 X 和 >》 方向 生成 。 生 成 的 高 斯 核 涉 及 高 斯 函数 
的 卷 积 和 高 斯 函数 的 一 阶 导 数 。 


(a) (0) ' 
g(x)= e 2 (11) 
图 12 原始 Lena 图 像 和 高 斯 模糊 处 理 后 的 Lena 图 像 OV27 
网 高 斯 函数 的 一 些微 分 如 下 所 示 : 
通常 来 说 ,一 个 函数 三 ( x, y ) 就 张 量 积 而 言 , 可 以 表示 为 
6 1 -2x -i 
;3 h 大 一 e 20 (12) 
f(y)= g(x)h(y) () i 
这 里 不 用 应 用 2D 内 核 ， 采 用 可 分 离 的 滤波 方法 来 计算 一 
维 内 核 沿 区 和 方向 的 梯度 。 i 
高 斯 函数 是 可 分 的 ,可 以 分 解 为 两 个 1D 高 斯 函数 的 乘积 。 Oe (13) 


[人 [全 
f(x,y)= = py ©9) g'(x)=—x*g(x)*o? (14) 
co) 使 用 生成 的 内 核对 图 像 执行 高 斯 平滑 。 


高 斯 滤波 器 的 重要 性 质 是 满足 不 确定 关系 的 唯一 滤波 器 。 ep tin 
AxAw=>1/2 (10) 
其 中 ，Ax 和 Aw 分别 是 空间 和 频 域 的 方差 。 独 特 的 属性 分 别 Hx= g(x)*g'(x) 1 
在 空间 与 频 域 之 间 的 本 地 化 冲突 目标 之 间作 出 最 佳 的 折 中 。 对 Hy= Hx ~ 
于 过 滤器 ， 高 斯 核 的 响应 在 无 限 域 和 大 部 分 领域 都 是 非 零 的 ; 


它 是 非常 小 的 ， 因 为 是 指数 形式 。 
2.4 改进 高 斯 一 拉 普 拉 斯 算 子 的 边缘 检测 

本 文 提出 的 改进 高 斯 一 拉 普 拉 斯 梯度 边缘 检测 步骤 总 结 如 
有 


可 和 输入 图 像 可 以 是 灰 度 或 彩色 图 像 。 变 量 sigma 用 于 确定 
两 个 方向 上 的 高 斯 核 。 在 高 斯 梯度 边缘 检测 方法 中 选择 适当 的 
sigma 值 。 图 13 为 输入 的 Lena 图 像 ， 图 14 为 灰 度 处 理 的 图 像 
以 及 高 斯 滤波 处 理 的 图 像 。 


(a) (b) 


pa 


15 ”原始 Lena 图 像 在 YX 和 了》 方向 上 的 高 斯 核 响应 
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qd) 对 生成 的 高 斯 平滑 进行 边缘 检测 。 


Edge Output 


到 的 输出 具有 分 割 效果 ; 对 于 Go 二 2， 发 生 分 割 ， 最 为 合适 的 
效 的 结果 。 边 缘 检 测 获 得 的 


O 二 1.5 时 , 产生 最 有 
所 示 。 


现 


16 原始 Lena 图 像 


其 中 : Gx 和 Gy 征 局 站 


(b) 


= abs(Gx)+abs(Gy) 
平滑 版 本 的 图 像 。 对 了 


和 改进 高 斯 一 拉 普 拉 斯 算 子 边缘 检测 处 


理 后 的 Lena 图 像 


3 ”实验 结果 


3.1 实验 设置 


在 本 节 中 ， 为 了 评估 所 提出 的 基于 改进 的 高 
高 斯 梯度 算法 的 性 能 , 在 Imagenet 数据 集 [15] 中 选取 了 10 幅 
像 进行 了 实验 ， 将 本 文 提出 的 高 斯 梯度 边缘 检 洲 


缘 检 测 器 进行 比较 。 图 


其 中 ， 图 17 (a) 为 原 图 
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图 17 显示 了 常规 边缘 检测 器 和 本 文 所 提出 的 高 斯 梯度 边 
缘 检 测 器 输出 结果 。 所 提出 算法 的 性 能 由 Imagenet 数据 集 图 像 
的 峰 品 比 (peak signal to noise ratio，PSNR)06~18 和 原 图 像 与 
处 理 图 像 之 间 均 方 误差 (mean square error，MSE)08-20 等 性 能 
度量 指标 进行 评估 。PSNR 是 最 普遍 、 最 广泛 使 用 的 评 鉴 画 质 
的 客观 


测 法 ， 常 用 做 图 像 压缩 等 领域 中 信号 重建 质量 的 测量 
方法 。 
是 原 图 


Er 


注 叶 地 六 认 全 


常 简单 地 通过 均 方差 (MSE) 进行 定义 。 其 中 ,MSE 
与 处 理 图 像 之 间 均 方 误差 。 


RR? 
PSNR =101o -一 -一 16 
a | (16) 


MSE = Dnll, (M,N)-1(M, NT 
M,N 

其 中 :J (M,N ) 代表 边缘 检测 输出 ，7o (M,N ) 代表 地 面 
真实 图 像 。 
3.2 ”实验 结果 

将 几 种 边缘 检测 算法 与 提出 的 方法 在 Imagenet 数据 集中 
选取 了 10 幅 图 像 进 行 了 实验 ,并 对 结果 进行 了 归纳 统计 。 表 1、 
~3 分 别 归纳 了 几 种 边缘 检测 算法 的 PSNR 值 、MSE 值 以 及 计 
算 时 间 ， 并 相应 给 出 了 折线 图 进行 直观 比较 。 几 种 边缘 检测 算 
法 的 PSNR 值 对 比如 图 18 所 示 ， 几 种 边缘 检测 算法 的 MSE 值 
对 比如 图 19 所 示 。 图 18、19 中 的 性 能 度量 图 显示 ， 提 出 的 改 


(17) 


沪 男 染 


17 为 其 中 一 个 图 像 上 的 


像 ; (b) 为 Sobel 算 子 检 涡 


Roberts 算 子 检测 的 结果 
高 斯 一 拉 普 拉 斯 算 子 检测 的 结果 ;， (9 为 改进 高 
子 检测 的 结果 。 该 仿真 ! 


; (d) 为 Prewitt 算 子 检 洲 


10 操作 系统 和 64 位 处 理 器 上 运行 。 


(f) 


图 17 其 中 一 个 图 像 上 的 边缘 检测 


MATLAB 2012a 软件 完成 。 
关 参 数 如 下 : Intel Corei3 @3.30 GHz、4GBRAM, 在 Windows 


疆 日 


进 高 斯 一 拉 普 拉 斯 方法 具有 高 PSNR 和 低 MSE。 算法 的 计算 时 
间 如 图 20 所 示 。 提 出 的 改进 高 斯 一 拉 普 拉 斯 方法 的 计算 时 间 
各 高 于 其 他 方法 ， 但 是 在 边缘 跟踪 的 质量 方面 ， 提 出 的 方法 优 
于 其 他 常规 技术 。 其 中 Sobel、Roberts 和 Prewitt 算 子 这 三 种 梯 
度 算 子 属 于 一 阶 微分 算 子 ， 对 于 噪声 及 光照 不 均匀 等 情况 的 存 
在 ， 不 能 完整 地 勾画 出 物体 的 边缘 。 而 高 斯 一 拉 普 拉 斯 在 平滑 
噪声 方面 的 效果 依然 不 理想 。 提 出 的 改进 高 斯 一 拉 普 拉 斯 方法 
则 很 好 地 解决 了 这 两 个 问题 。 对 于 PSNR 和 MSE 的 计算 ， 
Imagenet 数据 集中 提供 了 地 面 真 值 图 像 。 


表 1 几 种 边缘 检测 算法 的 PSNR 值 


高 斯 一 拉 普 改进 高 斯 


图 片 Sobel Roberts Prewitt 
拉 斯 拉 普 拉 斯 

图 片 1 17.34 17.45 17.06 17.62 17.98 
妈 片 2 17.28 17.35 17.03 17.55 17.87 
图 片 3 18.45 18.56 18.31 18.87 18.98 
到 片 4 19.04 19.13 18.85 19.29 19.34 
图 片 5 22.34 22.57 22.53 23.04 23.36 
到 片 6 19.34 19.86 19.21 19.88 20.42 
图 片 7 21.08 21.32 21.00 22.41 22.63 
到 片 8 17.89 18.08 17.45 18.33 18.56 
图 片 9 23.67 23.87 22.89 24.02 24.56 
图 片 10 20.73 20.86 19.91 20.85 21.03 
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1 五 
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图 y es EE 
Te 斯 梯度 算法 ,与 常用 的 梯度 算 子 相 比 ,取得 了 令 人 满意 的 结果 。 
图 上 请 4 六 全 壮 4 Sr 、， 村 、 > ; 
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氏 TFI ee 、 > 二 、 
0 斯 梯度 算法 的 性 能 ， 在 Imagenet 数据 集中 选取 了 10 幅 图 像 进 
医 | 2 、 本 
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人 人 人、 分 割 在 远程 医疗 应 用 中 对 分 析 感 兴趣 区 域 起 着 至 关 重 要 的 作用 
40 ee Cee a 
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